CONVERGENCE OF NONLOCAL THRESHOLD DYNAMICS 
APPROXIMATIONS TO FRONT PROPAGATION 



LUIS A. CAFFARELLI AND PANAGIOTIS E. SOUGANIDIS 

Abstract. In this note we prove that appropriately scaled threshold dynamics-type algorithms 
corresponding to the fractional Laplacian of order a € (0, 2) converge to moving fronts. When a ^ 1 
the resulting interface moves by weighted mean curvature, while for a < 1 the normal velocity is 
nonlocal of "fractional-type." The results easily extend to general nonlocal anisotropic threshold 
dynamics schemes. 



Introduction 

We study here the convergence of a class of threshold dynamics-type approximations to moving 
fronts. Although the arguments extend easily to general anisotropic kernels to keep the presentation 
simple, here we concentrate on a particular isotropic case, namely, fractional Laplacian of order 
a. € (0,2). The resulting interfaces move either by weighted mean curvature, if q E [1,2), or by a 
nonlocal fractional-type normal velocity, if a G (0, 1). 

Threshold dynamics is a general term used to describe approximations to motion of boundaries of 
open sets in by "measuring" interactions with the environment. The general scheme we consider 
here is described as follows: 

Let i7o be an open subset of with boundary Fq. The goal is to come up with an explicit 
approximation evolution scheme with time step /i > 0, so that, as nh — > t, the "approximate" front 
r^^, which is the boundary of an open set Q^nh identified as the level set of a sign-function, converges, 
in a suitable sense, to a moving front Ft, the boundary of an open set fit, and to identify the limiting 
velocity. For each ra G N, let 

where 

(0.1) n;,(.,0) = lno-Il«5 in , 

and, for n ^ 1, 

(0.2) uk{-, {n + l)h) = sign (A * Uh{-, nh)) . 

Here sign(t) = 1 if t > and —1 otherwise, 1a denotes the characteristic function of A C M^, and, 
for a € (0,2), 

(0.3) Jh{x) = paix,aa{h)) , 
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2 CafTarelli and Souganidis 

where pa is the fundamental solution of the fractional Laplacian pde 
(0.4) Wt - L'^W = in x (0, oo) , 

with 

(0.5) L"VF(x) = j \y-x\-^^+''^{W{y)-W{x))dy . 

Hence, at each time step, we solve the equation (|0.4p with initial datum 

W{-,0) =Uh{-,nh) in , 
and for time (Ta{h). Then we define Uh{--, {n + l)h) by 



Uh{x, (n + l)h) 



1 if W{x,ao,{h)) > 0, 
— 1 otherwise. 



The algorithm generates functions Uh{-,nh) and open sets defined by 

^nh = ■ Jh*Uh{-,{n-l)h){x) >0} and Uh{- , nh) = l^h -h^hy 

nh ^ nh' 



in 



We prove that, when /i — > 0, the discrete evolution Fq = dO.^^ converges, in a suitable 

sense, to the motion Fq Ft with nonlocal fractional normal velocity, if a G (0,1), and normal 
velocity equal to a multiple of the mean curvature, if a G [1, 2). 

If the kernel J is a Gaussian, it is a classical result that the algorithm generates movement by 
mean curvature — see below for an extensive discussion and references. 

To state the main result we recall that the geometric evolution of a front F^ = dVtt with normal 
velocity v{Dn,n,Qt), starting at Fq = d^o, is best described by "the level set" partial differential 
equation 



(0.6) 



ut + F{D'^u,Du,{u{-,t)^u{x,t)},{u{-,t)Su{x,t)} = in x (0, oo) , 
u = g on X {0} , 
with g such that 

(0.7) no = {xeR^ : g{x) > 0} and Fq = {x G : g{x) = 0} , 



and 



F(Z)^u, Du, {u{-, t) ^ u{x, t)], {n(-, t) ^ u{x, t)]) 
Du Du 

= \Du\v{-D{— -),--—, {ui-,t) ^ uix,t)},{u{-,t) < u{x,t)}) . 
\-lJu\ \Du\ 

The basic fact of the level set approach is that the sets Of = {x G : n(x, t) > 0} and 
Ff = {x G M.^ : u{x, t) = 0} are independent of the choice of the initial datum g provided the latter 
is positive in and zero in Fq. 

The weighted mean curvature motion corresponds to the level set pde 

(0.8) ut - Ca tr(/ -Du(g) Du)D^u = in x (0, oo) , 
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where, for p G \ {0}, p = p/\p\, while the equation corresponding to the nonlocal motion is 
(0.9) ut-Ca\Du\ / {l+{u{x + y,t) - u{x,t)) - r{u{x + y,t) - u{x,t)))\y\-^^+'''> dy = , 



where l"*" and 1~ denote respectively the characteristic functions of [0, cxo) and (— oo,0) and, in both 

cases, Ca is an explicit constant specified later in the paper. 

Although the heuristic meaning of (10. 8p is well known, some discussion about (10. 9p is in order 
It is implicit in (|0.9|) that there are sufficient cancellations in the term l"'"(n(x + y,t) — u{x,t)) — 
+ y, t) — n(x, t)) to compensate for the lack of integrability of the kernel y ^ at the 

origin. Indeed if we write the integral in polar coordinates. 



/■OO f- 



we see that the spherical integral measures the "deviation" of dQ, where Q = {y £ : u{y + x,t) ^ 
u{x, t)}, to be perfectly balanced between positive and negative parts, a deviation that infinitesimally 
is given by the mean curvature of d^. Observe also that if the surface {x € : u{y-\-x, t) = u{x, t)} 
is smooth and \Du\ 7^ on it, then an integration by parts leads to 

, Co /■ y Du{x + y,t) 

a J{ym^:u{^+y,t)=H^,t)} \y\ \Du{x + y,t)\ 

where dS is the (A^ — l)-surface measure. 

To state the main result of the paper we recall that, given a bounded sequence (uh{-,nh))n£N of 
bounded functions, the "half relaxed" limits u* and are defined by 



(0.10) 



u*{x,t) = \im sup* Uh{x,t) = \im sup Uh{y,nh) , 

y—>x , nh^t 

and 

U:^{x,t) = liminf^ u/i(x, t) = lim inf n/j(y, n/i) . 



It is immediate that ^ u* and, more importantly, if liminf^ Ufi = limsup* Uh, then, as /i ^ 0, 
Ufi ^ u locally uniformly. 
The result is: 

Theorem. Assume that Tq = dQo = (9(M^ \ Oq) and consider the family (u/j(-, n/i))„gi!} defined by 
(lOTD and dn^I) with aa{h) = if a e (1,2), h = o-2(/i)j In cJi(/i)| if a = I, and aa{h) = , 

if a e (0, 1). Let fit = {x G : u{x,t) > 0} and Tt = {x e R'^ : u{x,t) = 0}, where, for some 
uniformly continuous g such that ^0 = {x : g{x) > 0} and Tq = {x € : g{x) = 0}, u is the 

unique uniformly solution of (j0.6p with F given by (jO.Sp . if a^l, and (j0.9p . ifa<l. Then 

liminf^M/j = l in 0,t <ind lim sup* n/,, = — 1 in (Qf^'^tY ■ 

The Theorem asserts that the scheme characterizes the evolution of the front by assigning the 
values 1 inside the region {u > 0} and —1 outside the region {u < 0}. Whether the regions where 
Uh converges to 1 and —1 are exactly the regions inside and outside the front respectively depends 
on whether the front develops an interior or not. Interior are regions (patches) of positive measure 
where u = 0. Actually the answer is yes if and only if no interior develops. 

We have: 
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Corollary. // Ut>ort x {t} = d{{x,t) : u{x,t) > 0} = d{{x,t) : u{x,t) < 0}, then the set = 
U„gN(r^/, X {nh}) converges, as n ^ oo, to F = Uj>o(rt x {t}) in the Hausdorjf distance. 

The strategy of the proof of the Theorem is similar to the one in Barles and Georgelin [BGj . 
which is based on the general scheme developed by Barles and Souganidis |BS1] for convergence 
of monotone, stable and consistent approximation to viscosity solutions. Once the correct scaling 
<Tq-(/i) is identified, the main step is to prove the consistency of the scheme. The key difference 
between previous results and the one here is that all previous works considered kernels that were 
either exponentials or compactly supported, while in the case at hand they only have a prescribed 
power decay. 

Threshold dynamics schemes are used in probability and, in particular, percolation theory to 
find asymptotic shapes. We refer to Gravner and Griffeath [GGj and the references therein for a 
discussion and results from the probabilistic point of view. 

In the context of moving fronts, Bence, Merriman and Osher [MBO| introduced a scheme to com- 
pute mean curvature motion by iterating the heat equation. Evans [E] and Barles and Georgelin 
|BG] provided proofs for the BMO-algorithm. Some extensions to other isotropic kernels were con- 
sidered by Ishii [li], while Ishii, Peres and Souganidis jlPSj studied general anisotropic schemes with 
compactly supported kernels, and Slepcev [S] proved the convergence of a class of nonlocal threshold 
dynamics. Recently, Da Lio, Forcadel, and Monneau [DFM] studied the convergence, at large scales, 
of a nonlocal first-order equation to an anisotropic mean curvature motion. Although related to, the 
results of [DFMj are different than ours. Nonlocal operators and, in particular, fractional Laplacians 
of order a = 1 are often used to describe dislocation dynamics by line tension terms deriving from 
an energy associated to the dislocation line. We refer to Garroni and Miiller |GM1| . |GM2| for a 
variational analogue of what we are doing here for a = 1. For the case a < 1/2, the stationary 
solutions of (jO.gp must satisfy an integral "zero mean curvature" equation, that can be obtained 
as the Euler-Lagrange equation of a minimization process in the Hilbert space H"^"^ of functions 
with "a/2 derivatives in L^." A regularity theory of such surfaces, similar to the classical theory of 
"boundaries of sets with minimal perimeter," is being developed by Caffarelli, Roquejoffre and Savin 
|CRS| . Finally, Imbert and Souganidis |ImSj studied recently the onset of fronts at the asymptotic 
limit of fractional integral equations with reaction terms. 

Nonlocal phase transistion, models were proposed by Chen and Fife [CFj . Giacomini and Lebowitz 
[GLT] . pL2] . DeMasi, Orlandi, Presutti and Triolo [DOPTj . and DeMasi, Gobron and Presutti 
|DGP| . in the Landau-Ginzburg context of mean field theory for statistical mechanics. All the 
above references assume, however, fast enough decay or compact support for the diffusion kernels 
to guarantee an infinitesimal curvature condition for the limit. The connection between these non- 
local equations and the underlying stochastic Ising systems and moving fronts was established by 
Katsoulakis and Souganidis |KSH |KS2] and Barles and Souganidis |BS2j . 

The paper is organized as follows: In Section 1 we present some preliminaries. The proof of the 
Theorem begins in Section 2. The key step in the proof, i.e., the consistency of the scheme, is 
presented in Section 3. 

1. Preliminaries 

We recall some basic facts from the Crandall-Lions |CL| theory of viscosity solutions that we will 
be using in the paper. We begin with the definitions. Since the two equations (jO.Sp and (|0.9p are of 
different nature, i.e., local versus nonlocal, we give two separate definitions. 

Definition 1.1. An upper semicontinuous (use) (resp. lower semicontinuous (Isc)) function u is a 
viscosity suhsolution (resp. supersolution) of (jO.Sp if and only if, for all (/) € C^(M^ x (0, cxo)), if. 
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for any maximum (resp. minimum) point {x,t) € M.^ x (0, +00) of u — (j), 

(1.1) (f)tSCati{I -'D^(g)'D^)D'^4> if \D(t>\^0 or 4>t S if \D(t>\ = andD'^4> = 0, 
(resp. 

(1.2) (l)t'^CatiiI -'D4>®'D4>)D^(j) if \D(j)\ y^O or cpt if \D(j)\ = 0) andD^(j) = 0. 

To make precise statements for the nonlocal motion, it is necessary to introduce some additional 
notation. To this end, for u : x (0, 00) ^ M and A C M^, let 



(1.3) 



l+{v{y + x,t)-v{x,t)) - r {v{y + x,t) - v{x,t))\y\-^^+'^^ dy , 
l+iv{y + X, t) - v{x, t)) - ^iv{y + x, t) - vix, t))\y\-^^+''^ dy , 
l+(z;(y + x,t) - v{x,t)) - l~{v{y + x,t) - v{x,t))\y\~^^^°'^ dy 
t^{v{y + x,t) - v{x, t)) - lL(u(y + x,t) - v{x, t))\y\~'^^^'^^ dy 



I[v]{x,t) = 
L[v]{x,t) = 
lA[v]{x,t) 
l^[v]{x,t) 



where IL|_ and L denote respectively the characteristic functions of (0,cx)) and (— oo,0]. 
We rewrite the level set pde ()0.9p obtained for a < 1 as 



(1.4) 

We have: 



Ut - Ca\Du\I[u] =0 in x (0,oo) . 



Definition 1.2. A locally bounded use (resp. Isc) function u is a viscosity subsolution (resp. super- 
solution) of ()1.4p if and only if, for all (j) G C'^(M^ x (0, 00)), if for any maximum (resp. minimum) 



point {x, t) € M X (0, 00) of u — <f>, and for any ball 



C 



centered at (x, t), 



(1.5) 

(resp. 

(1.6) 



Few comments are in order here. Firstly, the difference between (II. Sp and (jl.6p is not a typo. It 
is actually necessary to guarantee the well posedness and, in particular, the stability of the solution 
— see [S] for a discussion of a similar problem. Secondly, it turns out (see Barles and Imbert (Bl] ) 
that Definition 1.2 is actually independent of 6. Hence, we may assume in the proofs that 6 is either 
fixed depending on (j), or, even, that 5 —> in an appropriate way. 

It is well known that the initial value problem (jO.Sp has a unique uniformly continuous solution — 
see, for example, Barles, Soner and Souganidis |iBSSj and Ishii and Souganidis [IS] for general results. 
The well-posedness of uniformly continuous solutions of the initial value problem for (jO.Op . which 
follows along the lines of the analogous result for integro-differential operators, has been studied 
recently by Imbert [Imj. 

It turns out, however, (see [BSSj for a general discussion) that very weak, e.g., discontinuous, 
viscosity solutions of (jO.Sp and (|0.9p may not be unique. The uniqueness is very much related to the 
issue of the development of interior. 
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We say that an evolving front (i7f,r(,]R^ \ Clt) does not develop interior if, for all t ^ 0, 

(1.7) U {Tt X {t}) = d{{x, t) : u{x, t) > 0} = d{{x, t) : u{x, t) < 0} . 
t>o 

There are several sufficient conditions on (f^o, Tq, \17o) that imply that there is no interior (see 
|BSSj ). A general necessary and sufficient condition, which is related to the uniqueness of solutions 
of ()0.6p . is given in the next proposition. For its proof we refer to |BS2| . 

Proposition 1.3. For an open subset Qq o/ R^, let {Qti^t,^^ \^t) be the level-set evolution of 
(r2o,ro,M^ \ CIq) with normal velocity —F. 

(i) The no-interior condition (11.71) holds if and only if it holds for t = and the initial value problem 
(j0.6p with initial datum uq = — IlKiv\no ^'^^ unique discontinuous viscosity solution. 

(ii) // ()1.7p fails, then every use subsolution (resp. Isc supers olution) w of (|0.6p with tf(-,0) ^ 
kio - %g (^^^P- ^(^ 0) = - ki-^J satisfies, in x (0, +oo), 

As far as the kernels pa are concerned we recall that, for (x, t) € x (0, oo), 

(1.8) p„(x,t)=t-^/-P„(xt-i/") , 

where, for some C]\f^a > and all x € M.^ , the kernel Pa{x) = Pa{x, 1) satisfies 

(1.9) ^ P^(x) ^ C7v,a(l + |a;r+°)~^ and \DPa{x)\ ^ CN,a\xf~^~^"{l + Ixf^''^^ ■ 
We will also use here that, both locally uniformly in \ {0} and in L'^{'K^), 

(1.10) lmit-Va(-,0 =Pa(-) , 

where, for some Cjy^a > 0, 

(1.11) Pa(x) = C7^,„|x|-(^+") . 

2. The proof of the convergence 

The main step of the proof of the Theorem is 

Proposition 2.1. The functions lim sup* Uh and lim inf^ tt/^ are, respectively, viscosity subsolutions 
and supersolutions of (10. 6p for the F specified in the Theorem. 

We postpone the proof of Proposition 12.11 and we proceed with the 

Proof of the Theorem. Let u be as in the statement and denote by sign* and sign^ the use and the 
Isc envelopes respectively of the sign function in R. 
The functions (see [BSSj for the proof) 

fl if n^O, fl if u>0, 

sign*(M) = < and sign^(M) = < 

[-1 if n<0, [-1 if u^O, 
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are respectively the maximal use subsolution and the minimal Isc supersolution of ()0.6p with initial 
datum lop — l^c — recall that = {x : g{x) > 0}. Therefore any subsolution (resp. supersolution) 
w of (j0.6p with the same initial datum satisfies 

w ^ sign* u (resp. w ^ sign^ u) in x (0, oo) . 

It then follows from Proposition 12.11 that 

(2.1) limsup*n/j ^ sign* n and liminf^n/j ^ sign^ u in x (0, oo) . 
Since Uh takes only the values ±1, (12. ip gives 

liminf^u/i = 1 in Qt = {u > 0} and limsup*Uft, = — 1 in {u < 0} . 
The proof of the first part of the Theorem is now complete. □ 

The Corollary follows exactly as in Section 5 of |BGj . 
We continue with the: 

Proof of Proposition \2.1\ We only present the argument for u = limsup* u^. The claim for liminf^ 
follows similarly. 

Let (/> be a smooth test function and assume that (xq, to) € ^ (0, oo) is a strict global maximum 
point of n — 0. To avoid any technical difficulties, we assume that \\m.\ni(^y g-^^^ (j){y, s) = +oo. 

If either u{xq, Iq) = —1 or (xq, ^o) belongs to the interior of the set {u = 1}, the facts that u is use 
and takes only the values ±1 yield that u = —1 in a neighborhood of (xo,to)- Therefore, 

(2.2) D,P{xqM) = , D^^{xoM) ^ and Mxo,to) = . 

Next we assume that {xQ,to) £ d{u = 1}. It is standard in the theory of viscosity solutions that, 
since (p grows at infinity and (xo,to) is a strict maximum, there exists a subsequence {xh,nhh) such 
that, as /i — > 0, 



(2.3) 



u*h{xh,nhh) - (l){xh,nhh) = max^N^^^{ul - </>) , 
ul{xh,nhh) 1 and {xh,nhh) (xo,to) , 



where li^ is the upper semicontinuous envelop of Uh- 

Here we used the upper semicontinuity of and the fact that 

lim sup* Uh = lim sup* . 

Once again Ufi taking only the values ±1 and u*^{xh,nhh) ^ 1, as /i ^ 0, imply that, for h 
sufficiently small, u'^{xh,nhh) = 1. Moreover, for all x G and n € N, 

(2.4) Uh{x, nh) ^ 1 — (j){xh,nh) + (j){x, nh) . 

Indeed, if u*j^{x,nh) = —1, the inequality is trivially true, while if u'^{x,nh) = 1, then 

(j){x,nh) - <j){xh,nhh) ^ , 

and, therefore, 

Uh{x,nh) ^ sign*{(p{x,nh) — (l){xh,nh)) . 
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Recall next that 

Uh{;nhh)= sign* {S{h)uh{;{nh-1)){-)) in , 

where 

S{h)v = Jh* V . 

Therefore, 

Uh{;nhh) ^ sign* {S{h)ul{;{nh - in , 

and, hence, 

ul{;nhh) S sign*(5(/i)n*(-,K " 1)^)) m ■ 

Let X = Xh- Since 

1 = u*f^{xh,nhh) = sign* S{h){ul{-, [nu - l)h){xh) , 

we have 

S S{h){ul{; K - mi^h) = S{h) Sign* <(•, K - 
The definition of S{h) and nj^ taking only the values ±1 yield 

^ / + x^,(n/j - l)h) - u*f^{xh,nhh)) 

(2.5) ^ 

- + x/i, (rife - l)h) - ul{xh,nhh)))paiy,(Ta{h)) dy . 

The proof will then be complete if we show that ()2.5p implies, that, at (xo,to)i if Dcf) = 0, then 
(l)t S 0, or, if \D(I)\ / 0, (/>t ^ i^*/*, {^^C-, t) ^ n(xo, to)}), with F given by if a £ [1, 2) and 

^ if a G (0,1). 

This is exactly the consistency of the scheme, which we investigate in the next section. 

□ 



3. The consistency 

Since the argument is technical, it is necessary to look at several different cases depending on 
the range of a and whether D</> vanishes or not. As discussed earlier we only check the subsolution 
property. 

We begin with the case a € (0, 1). To this end let 
(3.1) Ca=(2f Pa{0,y')dy') CN,a , 

where Cjy^a is the constant in (jl.lip . 

Proposition 3.1. Fix a G (0,1), set aa{h) = and assume that (j2.5p holds. Then, for any 

hall Bs C centered at (xo,to); (|l-5p holds at (xo,to) with Ca given by (|3.ip . 

Proof. Throughout the proof, to make the notation simpler, we drop the explicit dependence of pa 
and aa on a, we write a instead of o"(/i), th = nh and (j)hiy, s) = 4>{y + x^,, s) — (j){xh, th). 
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In view ()2.3p . for any A C M^, we have 

' {y G : tx* (y + Xh,th - h) ^ ul{xh, t/,)} nAC{yeR^ : Mv, th-h)^0}nA, 



(3.2) 



and 

pN 



^ {yeR^ : uliy + Xh,th - h) < ul{xh,th)} n A ^ {y e : c^h{y,th - h) < 0} H A . 



Let a = 0j(xo,to)- Since, as /i — > 0, (ptixh,th) — > </>t(xo,to)) for any 7 > and sufficiently small h, 
which may depend on 5, we have 

(3.3) (p{- + Xh,th - h) S + Xh,th) - {a - j)h in Bg , 
and, hence, 

(j)h{-,th - h) S (l)h{-,th) - {a- -f)h in Bs , 

where now the ball Bs is centered at the origin. 

We consider next two different cases depending on whether \D(j){xo,tQ)\ vanishes or not. 

If \D(l){xo,to)\ ^ 0, then choosing 5 sufficiently we may assume that the level sets of (p in Bg near 
0(xo,io) are graphs of smooth functions with uniformly bounded derivatives. 

If 

Ih = y" i^'^iuliy + Xh,th - h) -ul{xh,th)) - T {ul{y + Xh,th - h) - ul{xh,th))p{y,a) dy , 
in view of (|3.2p . we have 

(3.4) IhSt + K, 
where 

4 = / i^'^iKiy + Xh,th - h) -ul{xh,th)) - l~{ul{y + Xh,th - h) - ul{xh,th))p{y,cr) dy 

and 

4= /" {l+i(l)hiy,th-h))-r{My,th-h))piy,a)dy . 
Jbs 

Since, as /i ^ 0, u*^{xh,th) u*{xo,to), it follows from (II. Sp . (ll.lOp and the upper semicontinuity 
of that 

(3.5) limcr^^I^ ^ CN,alB^Au\ixo,to) . 

h-*0 

Next we concentrate on the limiting behavior, as /i ^ 0, of for 6 sufficiently small. 
Using (j3.3p we find that 



I^^III^ + <'+ + IV^'- 

where 

Mi = / {i+iMy,th))-i'iMy,th)))piy,'j)dy 
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and 

]Vj'± = ± / {l^{My,th)-{a-j)h)-l^{cl)hiy,th)))piy,a)dy. 
JBs 

We treat each one of these terms separately. As far as is concerned, we observe that the 
Dominated Convergence Theorem and (jl.lOp yield 

lE^ a-^ < CN,a / (l+(<A(y + xo, to) - 0(xo, to)) 
(3-6) -r{<p{y + xo,to)-<t>{xo,to)))\y\-^''^''Uy 

= CAr,J^J(/)](xo,to) , 

provided we can show that the integrand in EI^ divided by is integrable. 
To this end, using radial coordinates and (II. 9p . we find 



Jo 



where 



Di{r)= I {l+{q){y + Xh,th)-Hxh,th))-l i(l){y + Xh,th)-(^{xh,th)))da . 

JdB(r) 

In general we know that Df^{r) is bounded for all r. However, it is a calculus exercise to check 
that the regularity of (j) and the fact that \D(j)\ ^ in i?^ — recall that 6 is sufficiently small and 
\D(p{xo,to) \ 7^ 0, yield, for some constant C > depending on 5, 



Hence, since a G (0, 1), 



Diir) £ Cr . 



no 
JO 



We continue now with the analysis of ±IV^'^. We show that, for a = 

(3.7) m±a-^w'^ S -ia-j)\D<Pixo,to)r\2 [ P{0,y')dy') . 

Combining (|3.3p and (|3.7p yields, after letting 7^0, (jl.Sp at least when \D(l){xo,to)\ ^ 0. 
We show ()3.7p when a > 0. The case a ^ follows similarly. To this end, given that \D(f){xo,to)\ 7^ 
0, by taking 5 even smaller if necessary, we may assume that (j)hi-,th) has the form 

(3.8) (l)h(y, th) = (3yy + {Any, y) , 

with Ah = 2D'^(j){xh,th) and Ph = \D(ph{xh,th)\ \D(l){xo,to)\, as h ^ 0. 
We then have 

{y G : ^ th)} = {y G : ^ yi + (Ay, y)} 
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and 

{y G : (a - j)h S M-.th)} = {y G : a^h ^ Vi + {Any, y)} , 

where 

= Ph^Ah and ah = Ph^{a - 7) . 

Consider the integrals 



Jr = ± / t^{(t)h{y,th)-{a--i)h)-l^{(t)h{y,th))p{y,a)dy . 

In view of (II. 8p and (II. 9p we have, for some q. > 

0^p(y,<T) ^<jC^^Jyj-(^+") in . 
Moreover, as /i — > and almost everywhere in y, 

\t^{(t)h{y,th)-{a--i)h)-i^{(i)h{yM)\^^ ■ 

Therefore 

lim = . 

Adding and subtracting J^'^ to (j3.4p we see that to conclude the proof we need to show that 



(3.9) 



Let 



lim ±a-i / {t^{<Ph{yM) - {a-l)h) - l^{<l)hiy,th)))piy,(T)dy 

a — >0 



< 



-(a-7)|I?(/>(xo,to)rW PM)dy' . 



V'(ct) = / {t+{(t)h{y,th)-{a--i)h)-l+{(t)h{y,th))p{y,'y)dy 



and 

^h{y) = {Ahy,y) . 

Since a = ha^^^" and 
we have 

and 

l+((/.,(aV"y,t^)) = l+(yi +al/"xl/,(y)) . 
It is therefore immediate that ■0(0) = 0. Hence to prove (jS.Op we need to find V''(0)- 
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In view of the above simplifications we rewrite tp as 

V'(ct) = / (l+(yi + a'/''^h{y) - ana) - + a^'''^ h{y)))P{y) dy . 

We state below the main step of the proof of (|3.9|) as a separate lemma. To this end, let 
/ : [0, oo) ^ M be defined by 

(3.10) f{a) = [ + F{y, a))P{y) dy . 

We have: 

Lemma 3.2. Let f be given by ()3.10p and assume that F € C^(]R^ x [0, oo)) and 
dl^.F e C7(M^ X [0,oo)). Then 



f{a)-f{0) = - / l+{yi+Fiy,a))dy,{dpF{y,p)P{y))dy 
Jo Jr^ 

(3.11) -/ l+iyi + F{y,a))dy,F{y,a)P{y)dy 

+ £ I r+{y^ + F{y,p))d,{dy,F{y,p)P{y))dy . 

We continue the ongoing proof and return to the proof of the lemma later. We use Lemma [X 
with 

F(y, a) = a^/'^iAhy, y) - a^o and F(y, a) = ^^/"(^y, y) , 
both of which satisfy the assumption of Lemma 13.21 In either case we find 

N 

dy,F{y,a) = 2ai/"(iw,iyi + • 

Therefore, 

lima-'[£ I l+{yi + F{y,p))dp{dy,F{y,p))P{y) - J 1+ {yi + F{y,a))dy, F{y,a)P{y) dy] = . 
Since a £ (0, 1), it is hence immediate that 

^^'(0) = ah [ l+(yi)9,P(y) dy = -aj, [ d{y,)P{y) dy = - , " ~ / P{y') dy , 

Jrn J^n \D(p{Xh,th)\ JrN-i 

where S denotes the usual Dirac mass. 

Next assume that \Dip{xo, io)| = 0. Since, as h ^ 0, Ph = {Dfl^ixh, th)\ 0, the previous argument 
does not work and it is necessary to look at several cases. It is clear from the definition of the solution 
that we only have to show that 

a = (j)t{xo,to) ^ . 

First we assume that, along some sequence /i — > 0, /?h 7^ and o-p^^ 0. In this case it is 
possible to repeat the previous argument. The main difference is that, instead of a fixed 5, here we 
use 5h = PhC with C = 2\\D^M. 
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We look at each of the estimates/hmits of the first part. We begin with I^. Using that 
we find 

For the Umit involving EI^ we observe that the integration takes place over a set such that 
DhC{ye Bs, : Phlvil S C{\yi\^ + |yf )} C {y G Bs, : Phlvil S C SM + C\yf} 
^{y(^Bs, : |yi| ^/3h(2C)-i|yf } , 



Similarly, it follows that 



limsupcj ^Ph\Jh^'^\ = 

h~*0 



The only thing left to check now is that 
h^o J 

PM)dy' . 



< 



-(a -7) 



This, however, follows exactly as before. Notice that, since a<l, cr-^/"/3^"^=crQ -^cr/?^^— >0, as 

The next case is that, along a subsequence h ^ 0, either /3/i = or cr/?^^ oo. Here we argue by 
contradiction and assume that a > 0, and, hence, a — 7 > for 7 sufficiently small. 
Arguing as in the beginning of the proof of the \D(l){xo,to) \ ^ case, we find 



(3.12) 



S CaS''' + (l+(/?/,yi + (Ay,y)-(a-7)) - t' {l3hyi + {Ahy^y)-{a-^)))v{y.<y) dy 

JBs 



, -1 (cj ^l5hy\^(y^~^{Ahy,y)-{a--i)))P{y)dy . 
It is clear that, as /i — > and a.e. in yi, 

^{a-^i3hy\^(y^'^{Ahy,y) - (a-7)) ^ and r (a^Vftyi + ^""H^fey, y) - (a-7)) ^ 1 • 

Hence dividing (j3.12p by g and letting ^ we get a contradiction to a > 0. 

The last case to consider is that, along a sequence /i ^ 0, 7^ and c/?^"*^ — > t with -£ > 0. 

Again we rewrite (|3.12p as 



(3.13) 



0^Ca5-+ / (l+(yi + ai/°/3-i(A^y,y)-(a-7)cT/?^-') 

- l-(yi + ai/"/3^-i(Ay, y) - (a - ^)a(il^))P{y) dy 
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Since cr^/"/?^^ = ^ 0, as h ^ 0, letting /i ^ and 7 ^ in (IXTHI) gives 

^ / (l+(yi - a£) - r (yi - a£))P(y) , 

which impUes, in view of the symmetry of P, that we must have a ^ 0. □ 

We continue now with the 

Proof of Lemma \3.S\ To keep the ideas clear we present a formal proof using Dirac masses etc.. 
Everything can, of course, be made rigorous considering smooth approximations of 1"'' and passing 
to the limit. We leave it up to the reader to do so. 
For /> > we have 

f'ip) = [ S{yi + F{y, p))dpF{y, p)P{y) dy 

dy, + F{yp)))dpF{y, p)P{y) dy - f 5{yi + p))dpF{y, p)dy,F{y, p)P{y) dy 

l+{yi+F{y,p))dy,[dpF{y,p)P{y)]dy- [ {dpl+ {yi + F{y, p)))dy,F{y, p)P{y) dy . 

Therefore 

fia) - /(O) = -£ I + Hy, P))dy, [dyF{y, p)P{y)] dy 

t+{yi + F{y,a))dy,F{y,a)P{y) dy 

+ r [ l+{yi+F{y,p))dp{dy,F{y,p)P{y))dy . 
Jo JR^ 

□ 

We proceed now with the case a € [1,2) where, with the appropriate choice of iTq,, we obtain a 
front moving by a weighted mean curvature. We treat two separate cases namely a € (1,2) and 
a = 1. The reason is that the former is more or less straightforward, while the latter requires a bit 
more delicate analysis due to the border-line integrability properties of the kernel pa- 

Let 



(3.14) 



We have: 



PaiO,y')dy' 



-1 

,2 



yiPc,{0,y')dy . 



Proposition 3.3. Assume a € (1,2) and set (Ta{h) = // (12.51) holds, then, at (xo,io)? we have 

(fTT]) with Ca given by (IXTiD if \D(j)\ ^ or (pt S if Dcj) = and D^cp = 0. 

Proof. As in Proposition 13.11 to simplify the notation we write a for aa{h), t^ for nh, and (phiy^t) 
for (l){xh + y,t) - (l){xh,th). 
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The starting point is again ()2.5p which, using ()2.3p and after a rescahng, imphes 
(3.15) S I (l+(</>/,(cTi/"y, th - h)) - l'{<Ph{cy^l''y, tn - h)))P{y) dy . 



Expanding (j)h, using that <j)h{^,th) = 0, (j)h & C^'^ and a = /i"/^, we find 



<Ph 

where 



Ph = D(/)h{xh,th) , Ah = ^D'^4>h{xh,th) and at = 4>ht{xh,th) ■ 

After a rotation and a change of variables, we may assume that ph = Ph^i with = \ph\. We 
denote by Ah the matrix we obtain from Ah after the rotation. The integration in (|3.15p is taking 
place over the sets Ch and C^, where 

Ch = {yeR'' : a^l'^duvx + y)-ah^ 0[[a^l^\y\^ + a''l^)[\yf + a^/")) ^ 0} . 

We argue now as in the proof of Proposition 3.1 and [BGj. The difference with the former is that 
a € (1,2) gives integrability at the origin. The difference with the latter is that here the kernel only 
has algebraic decay while in [BGj it is an exponential. 

Here we only present the argument if |L)(/)(xo, io)| 7^ 0. If |Z)0(xo,to)| = and -D^0(xo,to) = 0; 
it is necessary to look again at different cases as in the proof of Proposition 3.1 and [BGj . The 
argument is considerably simpler than the one presented in Proposition 3.1 since we do not need to 
consider special balls, etc.. We leave the details to the reader. 

Since, as /i — > 0, /3ft — > |-Di;/)(xo, to)| / 0, the /?ft's are strictly positive for sufficiently small h. 

Therefore 

Cft = {y G : yi + y, y) - + 0{a^l^ \y\ + a''l'^)i\yf + ci^/-)) ^ 0} . 

Finally using that, as /i — > 0, Ah — > A and Uh ^ a, where A is the rotated matrix A and 
a = (j)t{xo,tQ), we find 

= {y G : yi + a'/''Ph{{Ay, y) - a + 0{a'/''\y\ + a^/''){\yf + a^/") + o(l)(|y|2 + 1) ^ 0} . 
After all the above reductions we are left with the inequality 

(3.16) S J (l+(*h(y)) - r (M/ft(y)))P(y) dy , ' 

where, for y G M^, 

^hiy) =yi + aV"(/5^)-i[(iy, y)-a + 0(c7V"|y| + a2/")(|y|2 + a^/") + o(l)(|y|2 + 1)] . 

_0_ 

Let R = a q for some > 0. Then 

i^{^h{y))P{y)dy= [ i^{^h{y))P{y)dy+ f i^{^h{y))P{y)dy . 

JBr JB'i, 
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Using (jl.gp . we find, for some cm, a > 0, that 



/ t^{^h{y))P{y) dy ^ Cn,^ [ (1 + lyl')"^ dy S CN,aR-" 
Jm\BR J\y\^R 

Fix 7 > 0. For h small we liave ^ ^'^ in Br, where 



Hence 



We summarize the above, using that 

r^{^h{y))P{y)dyS J r^{^''{y))P{y)dy and J l{^''{y))P{y))dyS J l-{^h{y))P{y) dy 
in the inequality 

(3.17) S J (i+(yi + pPh'i'iy)) - + Pl3h'^{y)))P{y) dy , 

where p = a^^" and 

iPiy) = (i + 7ld) - a . 

Let 

f{p) = [ - r{yi+pPz'ij{y)))P{y)dy . 



The properties of P yield that /(O) = 0, therefore, as in the proof of Proposition 13. H we use 
Lemma to find /'(O). 

It is a straightforward computation to see that it yields to the inequality 

a S C^[tr(i - - 2-\N + 1)7] . 

An elementary linear algebra calculation yields 

tr(i - ii,i) = i tr(/ - D^OroT^o) «) D(f(^to))D^^{xo,to) , 

and, hence, after letting 7 — > 0, to the desired inequality. □ 

We continue with the case a = 1. The argument is very similar to the one above. There is, 
however, a technical complication due to the logarithmic integrability of the kernel pi. To deal with 
this difficulty, it is necessary to choose ai differently. 

Let 

(3.18) Ci = (2 / P^{0,y')dy')-^UN-i Hm (i?^+iPi(0, i?)) , 

where lon-i is the area of the unit sphere in M^~^ and, as always. Pi is defined by (jl.Sp . 
We have: 
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Proposition 3.4. Assume a = 1 and choose a so that h = af{h)\lnai{h)\. If ()2.5p holds, then, at 
{xo,to), dni) holds with Ci given by iKW if {Dcpl or (pt S 0, if Dc/) = and D'^cp = 0. 

Proof. We only discuss the case \Dcj}{xQ,tQ)\ / 0. The argument when \D(I){xq, to| = and D'^cj}{xQ,tQ) = 
0, is similar to the one in Propositions I3.l l and 13.31 hence we omit the details. The proof follows very 
closely the one presented for the case a € (1,2), hence, again we only sketch the main steps. The 
main difference/difficulty is the logarithmic integrability of the kernel. To circumvent this potential 
problem it is necessary to make a "more involved" scaling. 

To this end for each 6 > 0, arguing as before, we reach the inequality 

(3.19) 0^ / [l+{yi + Fiy,p))-r{yi + F{y,p))]P{y)dy + CN,i6-^(T 

where, for p = a\lna\, 

F(.y, P) = <y{{Ah + 7h Id)y, y) - pan 

with Ah = f^h^A, = Ph^7 and ah = Ph^l- 
Define 

f{p)= I {t+{yi + F{y,p)-t-{yi + F{y,p)))P{y)dy . 

Using the decay properties of P and the dominated convergence theorem we find easily that 
/(O) = 0. 

As before we need to calculate /'(O). Arguing formally — the calculation can be justified rigorously 
using regularizations of l"*" and 1~, etc., — we get 

fip) = 2 / 6{yi + F{y, p))dp{F{y, p))P{y) dy 

+ / {l+{yi + F{y,p))-l-{y^+F{y,p)))P{y)dT.''{5a-^)' 

where dTi^ is the surface measure on dB^^-i and / denotes differentiation with respect to p. 
Then 

/(p)-/(0) = Ip + Ip, 

with 

1, = 2 r [ S{yi + F{y, X))dx{F{y, X))P{y) dy dX , 

Jo J\y\£5a-^{\) 

and 

Hp = r / [{l+i5a-\X)yi + Fi6a-\X)y, A)) 

Jo J\y\ = l 

- r {6a-\X)yi + F{6a-\X)y, X))]P{6a-\X)y){6a-^)^-^d^ dX 
where dS is the surface measure on dBi and A = o"(A)| ln(T(A)|. 
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The growth of P and the fact that a'{X) = (| Inaj — 1)^^ yield, for some C > 0, the estimate 



and, hence. 



hm p"^|E„| = . 



Next we analyze Ip. We begin with the observation that 

dxF{y,X) = {Ahy,yy-a 

Hence 

T — -I- 



with 



and 

rp 



ll = -2a r [ 5{yi + F{y, \))P{y) dy dX 

Jo J||/|^o-A-i 



ll= /"(llnal-l)-! / S{yi+Fiy,X)){Ahy,y)P{y)dydX 

Jo J\y\S5a-'^ 



It is immediate that 



\iinp-H\ = -2a I PiO,y')dy' 

p^O Jrn-1 



while 



p-'ll= r\n{5a-'){\\na\-l)-\ln5a-')-' [ S{yi + F{y, X)){Af,y,y)P{y) dy dX 

Jo J\y\SSa-^ 

and, hence, 

limp-il2 = -J—- hm (-i- [ \yfP{0,y')dy)tv{A-A,,,) . 
p^o ^ [1\ — i) R^oo Inn J\y\^R 

Returning now to (j3.19p we find 
(3.20) OSp-' [ {l+{yi + Fiy,p))-riyi + Fiy,p)))Piy)dy + CN,ijr^^ . 

Letting p — > yields 

2a I P(0, y') dy ^ —1— lim -1- / |yf P(0, y') dy 
Jrn-i [N - 1) R-^oo In R J\y\<R 



(3.21) 



\y\SR 

^ LON-i lim (fl^+^Pi(0,i?)) . 



iV - 1 R^oo 

We may now conclude as in Proposition 3.2. □ 



Convergence of nonlocal threshold dynamics approximations to front propagation 



19 



References 

[BG] G. Barles and C. Georgelin, A simple proof of convergence for an approximation scheme for computing motions 
by mean curvature, SIAM J. Numor. Anal. 32 (1995), 484-500. 

[BI] G. Barles and C. Imbert, Second-order elliptic integro-difjerential equations: viscosity solutions, theory revis- 
ited, Ann. IHP (2007), to appear. 

[BSS] G. Barles, H.M. Soner and P.E. Souganidis, Front propagation and phase-field theory, SIAM J. Control Optim. 
31 (1993), 439 469. 

[BSl] G. Barles and P.E. Souganidis, Convergence of approximation schemes for fully nonlinear second order equa- 
tions, Asymptotic Anal. 4 (1991), no.3, 271-283. 

[BS2] G. Barles and P.E. Souganidis, A new approach to front propagation problems: theory and applications, Arch. 
Rational Mech. Anal. 159 (2001), 109-135. 

[CRS] L.A. Caffarelli, J.M. Roquejofre and O. Savin, in preparation. 

[CF] C.-K. Chen and P.C. Fife, Nonlocal models of phase transitions in solids. Adv. Math. Sci. Appl. 10 (2000), 
no.2, 821-849. 

[CL] M.G. Crandall and P.-L. Lions, Viscosity solutions of Hamilton- Jacobi equations. Trans. Amer. Math. Soc. 

277 (1983), no.l, 1-42. 

[DFM] F. Da Lio, N. Forcadcl and R. Monncau, Convergence of a non-local etconal equation to anisotropic mean 
curvature motion. Applications to dislocation dynamics. Journal of European Math. Soc, to appear. 

[DGP] A. DeMasi, T. Gobron and E. Presutti, Traveling fronts in non-local evolution equations. Arch. Rat. Mech. 
Anal. 132 (1995), 143-205. 

[DOPT] A. DcMasi, E. Orlandi, E. Presutti and L. Triolo, Stability of the interface in a model of phase separation, 

Proc. R. Soc. Edinb. 124A (1994), 1013-1022. 
[E] L.C. Evans, Convergence of an algorithm for mean curvature motion, Indiana Univ. Math. J. 42 (1993), 

635-681. 

[GLl] G. Giacomini and J.L. Lebowitz, Exact macroscopic description of phase segregation in model allows with long 
range interactions, I. Macroscopic limits, J. Stat. Phys. 87 (1998), no. 1-2, 37-61. 

[GL2] G. Giacomini and J.L. Lebowitz, Phase segregation dynamics in particle systems with long range interactions, 
11. Interfce motion, SIAM J. Appl. Math. 58 (1988), no.6, 1707-1729. 

[GG] J. Gravner and D. Griffeath, Threshold growth dynamics. Trans. AMS 340 (1993), 837-869. 

[GMl] A. Garroni and S. Miiller, T-limit of a phase-field model for dislocations, SIAM J. Math. Anal. 36 (2005), 
1943 1964. 

[GM2] A. Garroni and S. Miiller, A variational model for dislocations in the line tension limit, Arch. Ration. Mech. 
Anal. 181 (2006), no.3, 535-578. 

[Im] C. Imbert, Fractional mean curvature flows, preprint. 
[ImS] C. Imbert and P.E. Souganidis, in preparation. 

[Is] H. Ishii, A generalization of the Bence, Merriman and Osher algorithm for motion by mean curvature, in Cur- 
vature flows and related topics (Levico, 1994), vol.5 of GAKUTO Internat. Ser. Math. Sci. Appl., Gakkotosho, 
Tokyo, 1995, 111-127. 

[IPS] H. Ishii, G.E. Peres, and P.E. Souganidis, Threshold dynamics type approximation schemes for propagating 

fronts, J. Math. Soc. Japan 51 (1999), 267-308. 
[IS] H. Ishii and P.E. Souganidis, Generalized motion of noncompact hypersufaces with velocity having arbitrary 

growth on the curvature tensor, Tohuku Math. J. 47 (1995), no.2, 227-250. 
[KSl] M.A. Katsoulakis and P.E. Souganidis, Generalized motion by mean curvature as a macroscopic limit of 

stochastic Ising models with long range interactions and Glauber dynamics. Comm. Math. Phys. 196 (1995), 

61-97. 

[KS2] M.A. Katsoulakis and P.E. Souganidis, Stochastic Ising models and anisotropic front propgation, J. Stat. 

Physics 87 (1997), no. 1/2, 63-89. 
[MBO] B. Merriman, J. Bence and S. Osher, Diffusion generated motion by mean curvature. Computational Crystal 

Growers Workshop (1992). 

[S] D. Slepcev, Approximation schemes for propagation of fronts with nonlocal velocities and Neumann boundary 
conditions. Nonlinear Anal. 52 (2003), no.l, 79-115. 



